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Motivated by the measurements reported by direct detection experiments, most notably DAMA, 
CDMS-II, CoGeNT and XenonlO/100, we study further the constraints that might be set on some 
light dark matter candidates, Mdm ~ few GeV, using the Fermi-LAT data on the isotropic gamma- 
ray diffuse emission, fn particular, we consider a Dirac fermion singlet interacting through a new 
Z' gauge boson, and a scalar singlet S interacting through the Higgs portal. Both candidates are 
WIMP (Weakly Interacting Massive Particles), i.e. they have an annihilation cross-section in the 
pbarn range. Also they may both have a spin-independent elastic cross section on nucleons in the 
range required by direct detection experiments. Although being generic WIMP candidates, because 
,«— n they have different interactions with Standard Model particles, their phenomenology regarding the 

fvj isotropic diffuse gamma-ray emission is quite distinct, fn the case of the scalar singlet, the one-to- 

one correspondence between its annihilation cross-section and its spin-independent elastic scattering 
cross-section permits to express the constraints from the Fermi-LAT data in the direct detection 
exclusion plot, c° — Mdm- Depending on the astrophysics, we argue that it is possible to exclude the 
singlet scalar dark matter candidate at 95% confidence level. The constraints on the Dirac singlet 
interacting through a Z' are comparatively weaker. 

o 

^J I. INTRODUCTION 

Qh The cosmological data show evidence of dark matter (DM) [I], but the nature of the DM particle, if any, is still 

unknown. Recently, some interest has been taken in a light WIMP 1 , with Mdm ^S 15 GeV. This is largely motivated 
by the modulation observed by the DAMA/LIBRA and DAMA/Nal experiments [5] , the excess of events at low recoil 
energies measured by CoGeNT [3] , and an anomaly reported by the CRESST collaboration at recent conferences [3] . 
If interpreted in terms of Spin Independent (SI) elastic scattering of a DM particle, these experiments all point to a 
candidate with a mass in the few GeV range, and a SI cross-section er" = O(10 -40 ) cm 2 [5TllO|. Whether these results, 

^vj to which we will loosely refer as CoGeNT-DAMA, are due to dark matter or to a more mundane sort of background 

is a matter of debate (see e.g. [TTlll3| ). Furthermore other direct detection experiments have set exclusion limits that 
may exclude part or whole regions of the parameter space favoured by CoGeNT-DAMA (see for instance [7] and |14j). 
Most relevant for the low mass region, M DM < 15 GeV, are the CDMS-Si [15], CDMS-II [TB], XenonlO [HUTS] and 
XcnonlOO [TH] exclusion limits. However it is quite delicate to set robust exclusion limits in the low mass region, 
which in the direct detection experiments corresponds to small recoil energies 2 , and it may pay to look further for 

£ — . alternative ways to constrain a light WIMP. 

A most interesting possibility in the near future is to look for missing energy at colliders, and it has been shown 
that, depending on the assumed properties of the light WIMP, the limits by the Tevatron and the LHC may be 
competitive with direct searches for light candidates, Mdm ^ 10 GeV [25, 26] . The next best possibility is to consider 
^ indirect signatures that are specific to DM, i.e. with no or little background from known physics, like a monochromatic 
gamma-ray line on the sky, or a flux of energetic neutrinos from the Sun. The former signature depends very much 
on the dark matter candidate. The latter however is quite generic for a light WIMP relevant for CoGeNT-DAMA, at 
least because they must have a large scattering cross-section on matter, and thus may be efficiently captured in the 
Sun. Their subsequent annihilation at the centre of the Sun would release a flux of neutrinos which depends on the 
annihilation channels of the DM candidate, and may be constrained by Super-Kamionkande [SI 151 [2TH29] . 
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By a WIMP (Weakly Interacting Massive Particle) we mean a cold dark matter candidate, with an annihilation cross-section of the 
order of 1 pbarn, so that its relic abundance may be fixed by the standard thermal freeze-out mechanism. This potentially includes 
many candidates, whose interactions might have nothing (or little) to do with weak interactions. By a light WIMP, we mean here a 
candidate in the mass range between 1 and 15 GeV, substantially lighter than the usual suspects, like the neutralino. 

2 The exclusion limits may be very sensitive to parameters which are poorly constrained, see for instance the discussion revolving around 
the scintillation efficiency in LXe experiments l9l l20tf2H . 



The other indirect signatures, be it a continuous gamma-ray spectrum or extra antimatter in cosmic rays, are 
more prone to uncertainties, most notably because of large and not well-understood backgrounds from astrophysical 
sources. Nevertheless their measurement by the new generation of instruments (in particular Fermi- LAT) are expected 
to put very relevant constraints on a light WIMP. Indeed, although the annihilation of a light WIMP produces few 
particles (for instance, a 10 GeV candidate annihilating into bb quarks would give on average O(10) photons), this 
is more than compensated by the 1/M DM dependence of the fluxes, for fixed energy density pdm- This is of course 
well-appreciated, and, incidentally, the current analysis already give constraints on the annihilation cross-section that 
are relevant for generic light WIMP candidates (see for instance [5UH55] ). From the perspective of candidates suitable 
for CoGeNT-DAMA, this has been particularly emphasized in various recent works. To our knowledge, possible 
constraints from the isotropic diffuse gamma-ray emission or from the centre of our galaxy, are mentioned in [8] and, 
respectively, in [51136 1 157 ] . Also, constraints from Milky Way Dwarf Spheroidal galaxies (dSph) are discussed in (HUH]- 
In [pj in particular, the analysis of the Fermi-LAT collaboration of the gamma-rays flux from Milky Way dSphs [35J 
has been exploited to set a conservative constraint on a 10 GeV scalar singlet candidate, and to tentatively exclude 
lighter candidates, assuming a Navarro-Frcnk- White (NFW) [39, 40 profile for the dark matter in dSphs. Constraints 
from antimatter in cosmic rays, in particular anti-protons and anti-deuterons, have been discussed in |41H44j and 
more recently in [45] . where it is emphasized that candidates below 10 GeV are disfavoured by observations of the 
anti-proton cosmic ray flux. Another potentially important signature is synchrotron radiation from the electrons and 
positrons produced directly, or indirectly through inverse Compton scattering, in the annihilation of dark matter (see 
for instance [4"6T(4"9"] . and more in relation with CoGeNT-DAMA [50]). 

In the present article, we study further the constraints that may be set on light WIMP candidates related to 
CoGeNT-DAMA, based on the first-year Fermi-LAT data on the spectrum of isotropic diffuse gamma-ray emission [51] , 
Our motivation for doing so is manifold. The most obvious one is that we may anticipate that these observations are 
relevant for a generic light WIMP. In the analysis of [31] , the first-year Fermi-LAT data are used to set constraints 
on the annihilation cross-section of generic candidates, and this down to Mdm — 5 GeV. For the lightest candidates 
annihilating into t + t~, the constraints touch the region relevant for a WIMP, < av > an n ~ 3 • 10~ 26 cm 3 s _1 (see 
Figure 2 of Ref. [3T]). In the analysis by the Fermi-LAT collaboration [51], which covers generic candidates with 10 
GeV < Mdm < 1 TeV, the constraints appear to be more stringent and, depending on the astrophysical assumptions, 
some candidates may be excluded at 95 % confidence level (C.L.). These works also have a focus on more "popular" 
DM candidates, like the neutralino, Kaluza-Klein DM candidates or models designed to explain the possible excess of 
positrons in cosmic rays. 

The two analysis mentioned above essentially concur, and the differences in the strength of the constraints may be 
traced (albeit in a very convoluted way) to the use of distinct assumptions regarding the astrophysical uncertainties, 
in particular the distribution of inhomogeneities as a function of redshift. However, given the relevance of light 
WIMPs candidates for the anomalies reported by some direct detection experiments, and the current astrophysical 
uncertainties, we found of interest to make available yet another independent analysis, but with a focus on the low 
mass region, 5 GeV < M&M ^ 10 GeV. Moreover the data from Fermi-LAT have the potential to constrain candidates 
with a mass below 5 GeV, possibly down to Mom ~ 1 GeV. On one hand, going to lower masses is a tiny step, and, 
as a first approximation, one could simply extrapolate the existing constraints (from either |31j or [33]). On the other 
hand, the current direct detection experiments are not sensitive to such low mass candidates (be them unrelated to 
CoGeNT-DAMA) and, in our opinion, it is of interest to be as precise as possible. A further motivation of ours is that 
we would like to confront two specific, albeit quite generic, models of dark matter, in particular a real scalar singlet, 
interacting through the Higgs portal, and a Dirac fermion, interacting through a new Z' {i.e. the U(l)y portal). 
Although these models have similar pros (they are simple, minimal extensions of the Standard Model, hereafter SM, 
and they are directly relevant for CoGeNT-DAMA), and cons (SI elastic scattering with nuclei — which is challenged 
by Xenonl0/100 — and ad hoc mass scales and couplings), they have a distinct phenomenology with regard to isotropic 
diffuse gamma-ray emission. Specifically, these candidates have distinct couplings to SM light degrees of freedom (ie 
e + /e~~ and neutrinos) and thus have a distinct kinetic decoupling temperature, a parameter which is relevant to 
establish the mass of the smallest gravitationally bound DM objects that may form in the Universe, and consequently 
for the calculation of the diffuse isotropic gamma-ray flux. This implies that the scalar singlet (and its siblings with 
similar interactions) is sensibly more constrained than a candidate interacting through a Z', a feature which should 
be of interest for further model building. 

The plan of our article is consequently as follows. In section [11] we briefly expose our hypothesis regarding the 
calculation of the isotropic diffuse gamma-ray emission. As in [33] . we will suppose that the contribution from 
dark matter annihilation within the halo of our galaxy (and possible subhalos within) is negligible, and so, that the 
dominant contribution is due to the isotropic gamma-ray background radiation (IGRB). For the sake of comparison 
with other approaches, and to set our conventions regarding the astrophysical parameters, we discuss in section [TIT] the 
exclusion limits in general terms. Last but not least, in section [TV] we compare the outcome of our calculations for two 



minimal, toy models of dark matter which have been shown to be consistent with both CoGeNT-DAMA and WMAP, 
but have distinctive characteristics: a scalar candidate interacting through the Higgs portal, thus with '" Higgs-like" 
couplings to ordinary matter, and a Dirac fermion interacting through a Z', hence with "Z-like" couplings. Our main 
results arc given in Figures [6] and [Slfor the scalar singlet, and Figure [9] for the Dirac fermion. Finally we give our 
conclusion and prospects in sectionjv] 

II. THE ISOTROPIC EXTRAGALACTIC GAMMA-RAY FLUX 

In this section we summarize the formalism we have used to compute the extragalactic contribution to the gamma- 
ray from dark matter annihilation. We follow the approach of [52 and compare the result to other approaches in the 
next section. 

To begin with, we will work with the hypothesis that the contribution of dark matter annihilation within the halo 
of our galaxy to the diffuse isotropic gamma-ray emission is small compared with the extragalactic component, or 
IGRB, following the approach of !33 , to which we refer for future discussions on this point. 

The flux of gamma-ray today (in units of GeV _1 cm~ 2 s _1 sr _1 ) from DM annihilation at any redshift z is then given 
by [S3] 3 

rf$ 7 _ c <av > ann f°° , , 1 dN ~t K>2 ( j\ -t(e',o, z ') m 

dE 4vr 2M* M J H(z')(l + z') 4 dE> [ ' ' [) 

In this equation H(z)(l + z) = _ffo/i(z)(l + z) sets the relation between redshift interval and the proper distance 
interval, with h(z) — \/Oj^(l + z) 3 + fl\. For the cosmological parameters today we use the central value from 
the WMAP-7yrs measurements pQ: H = 70 .4 km/s/Mpc, fl M = 0.227, fl b = 0.0456 and fi A = 0.728. A factor 
of (1 + z)~ 3 accounts for the dilution of the number density of photons from expansion. Finally E' = E(l + z') 
is the energy of the photon at redshift z' for an energy E observed today. The factor dN 1 /dE is the spectrum of 
gamma-rays produced in the annihilation of a pair of DM particle, at a rate given by < av > ann , and depends on 
the particle physics candidate. The function B(z) takes into account the density profile of dark matter at a given 



redshift z and will be presented in Section II B Finally the parameter t(E' , z, z') is the optical depth (of a gamma-ray 
emitted at z' at energy E' and observed at z), which takes into account the absorption that may occur on the path of 
a gamma-ray For the case of a light WIMP (Mjjm ^,10 — 15 GeV), and focusing on the Fermi-LAT window, which 
has measured the isotropic flux between 200 MeV and 100 GeV, the Universe is essentially transparent. This is an 
interesting simplification as it removes one source of uncertainty on the calculated gamma-ray flux, so we discuss this 
aspect in some details in the following section. 

A. Optical depth 

For increasing energy, the dominant processes that may affect the propagation of a gamma-ray between redshift z 
and today are (see, for instance, |54j ) 

Photo - ionization 7 + H{He, He+) -> e" + H+(He+,He ++ ) , (2) 

Compton Scattering 7 + e~ — > 7 + e~ , (3) 

e + e~pair production on matter 7 + A — > A + e + + e~ , (4) 

Photon — photon scattering 7 + TekGrd ~~ ► 7 + 7 1 (5) 

e + e~pair production 7 + TBkGrd — ► e+ + e ~ • (6) 

To compute the processes (pi) and Q we have used the approach of [53], for Compton scattering we refer to |55j . 
to [56] for pair production on a background photon and, finally to [57) for the photon-photon scattering contribution. 
The results of our calculations are summarized in Figure [l] and they are in agreement with, for instance, [32j . There 
we show, in function of the observed energy today (dotted at z = 0), the redshift 1 + z' at which the optical depth was 
t(1 + z') = 1, such that the Universe may be considered to be optically thin (transparent) at lower redshifts. The red 
(gray) region corresponds to the energy window of Fermi-LAT, which shows that only the (Hj-pl processes are relevant 



3 The factor of 1/2 is for a self-conjugate DM candidate. In the case of a Dirac fermion which we will also consider, there is an extra 1/2 
factor 1531 . 



(see also [35]). The triangular region (in blue/dark gray in the plot) within the Fermi-LAT window gives the allowed 
redshifts at emission in function of energy for a candidate with mass Mom ^ 20 GeV. The fact that t(1 + z') « 1 
within the triangular means that the Universe is optically thin for these gamma-rays. However, we show in the plot 
only the effect from scattering or pair creation on the photons from the cosmoc microwave background (CMB), as is 
relevant for large redshifts, but, at lower redshifts, z % 6 — 10, we should also take into account the ultraviolet (UV) 
photons produced by the formation of the first stars (an effect not shown in Figure [I]). Unlike the CMB, there is 
much uncertainty regarding the distribution and spectrum of this extra background light. Here we follow |32j and 
also [55] and [5j5] , where the possible impact of the low z UV background is discussed. In |35j , three models for the UV 
background are considered. A glance at Figure 3 in Ref. [32] reveals that, for the case of a so-called realistic photon 
background, the effect is not relevant for a light WIMP. A similar conclusion may be drawn from the plots in Figure 
11 of [35], which shown that absorption due to UV photons is negligible within the range relevant for a light WIMP 
within the Fermi-LAT window. Absorption due to UV photons is more important according to the somewhat older 
model of [5_S], but, keeping in mind the fact that spectrum of photon produced by a, say, 20 GeV candidate typically 
peak at 1 GeV, we may safely conclude from this reference that absorption is essentially negligible, even in the most 
pessimistic scenario, if we concentrate on a light WIMP and the Fermi-LAT data. Thus, in the calculation of the 
gamma-ray flux, the only relevant astrophysical uncertainty is the distribution of dark matter halos as a function of 
redshift, which we now discuss. 

B. Enhancement factor from structure formation 

To compute the extragalactic contribution to the gamma-ray flux from dark matter annihilation, we need the 
distribution of dark matter halos as a function of redshift. There are different approaches to this complicated problem. 
Needless to say, this will be the main source of uncertainty on our results. Here we found convenient to follow a 
standard, semi-analytical approach based on the Press-Schechter formula. In particular we follow the notation of |52j . 
which is itself based on [60] . One should bare in mind that different approaches may give substantially different 
results. For a review on halo models and structure formation we refer to [611 162) . 

As the Universe expands, the initially small perturbations in the dark matter distribution evolve in a non-linear way 
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Figure 1: Contour plot of the optical depth with t(E' , z, z') = 1, where E' /(l + z') — E/(l + z). Here E is the energy observed 
at redshift z, so that E' is the emission energy at redshift z . The different lines correspond to the different observers. For 
each, they give the maximal possible redshift at emission: z = dotted line, z — 10 dashed curve, z = 100 long dashed line 
and z = 1000 solid line. The red (gray) region denotes the Fermi-LAT measured energy range. The triangular region gives the 
relevant (for observations today) redshifts for a DM candidate of mass Mdm = 20 GeV. 



into a complex network of structures made of dark matter halos of various scales. In the ACDM model for structure 
formation, dark matter halos are assumed to form hierarchically bottom-up, with small structures merging into larger 
halos, in a self-similar pattern. Under these hypotheses, the Press-Schechter [53] empirical formula gives an analytical 
expression for the abundance of dark matter halos of mass M at a given redshift z, on the basis of the predictions of 
linear perturbation theory for the power spectrum of perturbations: 

dn [2 p c n DM da(R, z) 5 C < 6* \ .. 

dM {z > M) = n^i — dM a2(jRi ,) ex Pv-2^(^)J ' (7) 

with p c being the critical density. The factor S c Eq. (17| is the critical over-density required in a spherical model for 
gravitational collapse. We take S c = 1.28 following |62| . The other key quantity is cr 2 (R, 0), which is the variance of 
the density field, obtained using linear perturbation theory, in a sphere of radius R containing a mean mass M. Its 
evolution with respect to redshift is simply given by the growth factor, using the standard factorization between scale 
dependence and redshift [51] . The variance of the density field on a scale R is related to the matter power spectrum 
P(k) through convolution with a spherical top-hat window function W(kR), with R 3 = 3M/4irp c £lM- To compute the 
linear matter power spectrum today, we have used CAMB [651 I66j , linearly extrapolated at very large wave- number 
k for the evaluation of Eq. Vty. Both P{k) and a are normalized by computing a in a sphere of R — 8/i _1 Mpc, with 
cr 8 =0.8 [531 EQ]. 

From the distribution of halo of mass M at redshift z, we may compute the flux using 

B 2 {z) = B 2 (z) + B 2 {z), 

/ 1 f dn fr vir (M,z) 

= ^ DM {l + zf(l+ ^ DAi{1 + z)3 JdM m (*,M)(l+«)'/ o dr4*rV(r,JO), (8) 

where the first term is the contribution from the smooth distribution of dark matter, which in practice is negligible, 
and p(r, M) is the dark matter density profile of dark matter in halo. The integral over p is calculated with a cut-off 
given by the virial radius r vir (M, z), which is usually expressed in terms of the so-called concentration parameter 
Cvir — r V i r (M, z)/r s where r s is the core radius which depends on the choice for the profile [57] • The virial radius of a 
halo of mass Mat a given redshift z is itself defined as the radius within which the mean density of the halo is A V i r 
times the smooth density p c ^M- In the ACDM model, this virial over-density is given by 

A , , 18^ 2 + 82(Q M (z) - 1) - 39(n M (z) ~ I) 2 , Q . 

which generalizes the standard result A V i r (z) = 18ir 2 ~ 200 of a CDM Universe. 

As is standard practice for primordial halos, we consider throughout the paper the NFW profile as a benchmark |39l 
140] . The NFW density profile has a simple analytic expression in terms of the concentration parameter 

f *-m = % i 1 - (TT^-f) ( Log(1 + ^ - (irS-) ) 2 ■ (10) 

Using the previous definitions, the integral over the halo mass function can finally be recast into the form 

A f 00 dn 

BW = -^ dMM — F NFW (c mr (z,M)), (11) 



where M m j n is the minimum halo mass (see Section III) so that Eq. (18j) becomes 



B(z)=p c n DM (l + z) 3 ^l + B(z). (12) 

This is the enhancement factor which appears in the extra-galactic flux, Eq. M. The dark matter distribution is 
described by the mean cosmological matter density B plus the sum over all average enhancements due to a halo of 
mass M, weighted over the mass function. The abundance and distribution of halos depend primarily on halo mass 
and the concentration parameter itself depends on the halo mass at a given redshift. Therefore, the range of M m [ n 
and the functional form of c V i r are the major sources of uncertainty in the determination of the diffuse gamma-ray 
flux and we will discuss our choices in Section [Ml 
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Figure 2: Example of E 2 dN 1 /dE gamma-ray spectra from the annihilation of a 10 GeV dark matter particle, as a function 
of the photon energy E. The solid blue line represents a pure annihilation into bb (BR b i = 100%), the green dashed curve is 
for a pure annihilation into t + t~ (BR t + t - = 100%), the red long dashed to a "Higgs-like" annihilation (BR T + T ~ ~ 9% and 
BR bb ~ 83%), while the black dot-dashed line refers to a"Z-like" annihilation (BR t ± ~ 3% and BR q q ~ 13%). 

C. Gamma-ray spectrum at production 

The last ingredient we will need to compute the contribution to the IGRB, is the spectrum of gamma-rays produced 
per DM annihilation, 



d/V 7 
~dE 



E 



dNl 

dE ! 



(13) 



where dN? /dE is the photon spectrum for the / annihilation final state with branching ratio BRf. The gamma-ray 
spectra from hadronization and final state radiation arc calculated with Pythia 8.1 [68 . In Figure [2] we show a few 
spectra representative of the particle models that we will analyze. The solid blue line is the final state radiation from 
annihilation into bb pairs, the green short-dashed curve comes from a BR = 100% into t + t~ . As is well known, in both 
cases, the main contribution comes from the decay of neutral pions into gammas, produced from the hadronization 
of the bb pair or from the semi-hadronic decay of the t + t~ . The most relevant feature for our analysis, which will 
make use of the lower energy edge of the Fermi-LAT data, is that the photon spectrum from t + t~ is harder than the 
one from bb quarks. The red long-dashed curve is the photon spectrum for light WIMP with "Higgs-like" coupling to 
SM degrees of freedom, namely with a branching ratio dictated by the Yukawa couplings. For instance, for a scalar 
dark matter particle with Mdm — 10 GeV, it corresponds to a BR bb = 83% and to a BR T + T - — 9%. The black 
dot-dashed line is the spectrum for a DM particle with "Z-like" couplings, in which case BR bb ~3x BRi± ~ 13%. 



III. GENERAL RESULTS REGARDING A LIGHT WIMP 



In this section we discuss, and compare to known results, our choice for astrophysical quantities like the concentration 
parameters and the minimum DM halo mass that may form. These parameters are defined in a non-unique way in 
the literature, and the choice of their parameterization will affect the limits on the annihilation cross-sections. Once 
the technical details are worked out, we present some generic constraints relevant for a light WIMP, with Mdm up 
to 20 GeV. 





Figure 3: On the left, boost factor 1 + B(z) as a function of the redshift z. On the right, zoom at small redshift values of the 
total enhancement factor B 2 (z)/h(z)(l + z) A . A NFW density profile for the dark matter is assumed. The thin green lines are 
given for the Cwmap concentration parameter |69] . while the thick black lines are for the Cpl concentration parameter [67] . 
The long dashed lines refer to M min = lCr s M , the solid to M min = 10" 6 Af Q and the short dashed to M min = 1O" 4 M . 



A. Uncertainties from structure formation 

Simulation of halos [39, 671 indicate that there is a strong correlation between the concentration parameter and the 
halo mass M. Among other properties, they are inversely related, e.g. higher mass halos are less concentrated. A 
model that fits the results of the simulations is: 



,{z,M) 



ir(0,M) 

1 + z 



(14) 



meaning that c V i r is inversely proportional to the redshift as the radius of a halo increases as the Universe expands. 
We use two different power law fits for c V i r . The first one (Cwmap hereafter) is the WMAP-5yrs best fit in [69] . 
which reads Log c vir = 0.971 — 0.094 Log [M vir / (10 12 h _1 M Q )]. The second one (Cpl from now on) follows a simple 
power-law behavior, c V i r oc M , as indicated by numerical simulations [67j and is normalized at the same mass 
scale as Cwmap- 



The other quantity that affects the final results, namely the upper bounds on < av >„„„, is the choice of the 



minimal halo mass, appearing as the lower limit in the integral of Eq. (111. It has been shown in [701 171] that free 



streaming and collisional damping of a WIMP with particles from the thermal bath lead to a small-scale cut-off in 
the DM density perturbation power spectrum. The size of the smallest gravitationally bound structure is set by the 
properties of the WIMP. After the DM particle abundance has frozen out in the early Universe, it may still stay in local 
thermal equilibrium by elastic scattering processes with SM particles, DM + f — > DM + /. As the Universe expands 
the elastic scattering rate decreases and the WIMPs kinetically decouple at some temperature, Tkd which depends on 
the WIMP properties. From this moment on the WIMPs can stream freely and wash out matter contrasts on small 
scales. For a candidate like a neutralino, the complete decoupling from the thermal bath typically happens when 
the temperature has dropped by a factor ~ 10 — 1000 with respect to the chemical decoupling temperature, namely 
Tkd ~ O(MeV) [7D1 [7T]. In the next section we will consider two distinct light WIMP candidates, a Dirac fermion with 
"Z-like" couplings to SM degrees of freedom, and a scalar candidate with "Higgs-like" couplings to ordinary matter, 
i.e. through Yukawa couplings. Here we estimate the temperature of kinetic coupling and the resulting estimates for 
the minimum halo mass, emphasizing the differences between the two models. 

In the case of a scalar DM interacting through the Standard Model Higgs, we could expect that kinetic decoupling 
occured at about the same temperature as chemical decoupling, typically Tkd ^ 100 MeV for the candidates of interest 
here. More specifically, for SM Yukawa couplings, the most relevant scattering is a priori with c, b and t~ and their 



antiparticles, which are non-relativistic, and thus have an exponentially suppressed abundance at the time of interest. 
As a first estimate for Tkd, we may balance the rate for elastic scattering with the expansion rate of the Universe. 
In this approximation kinetic decoupling is like chemical decoupling, but the density of dark matter is replaced by 
the density of target fermions in the collision rate. From this analogy, we estimate that kinetic decoupling occurs for 
2^kd = «v/21d = O(20) assuming a = 0(1O -36 cm 2 ), which gives Tkd ~ 250 MeV in the case of / = b/b, while for c/c 
and t + /t~ we have Tkd ~ 100 MeV. These temperatures are of the order of the chemical decoupling (fo) temperature, 
Tf ~ 250 — 500 MeV for Mjjm ^5 — 10 GeV. Things are however a bit more subtle, as momentum exchange typically 
requires that many collisions take place, N co \ > 1. From the analysis of Ref. [72], we have N co \ ~ Mdm/T > 1 or 
-V C oi ~ Mom I'm / ^ 1 for relativistic or non-relativistic targets respectively. Furthermore the interaction with /x + //i~ 
are small but still substantial, so we need a more precise formulation. For this, we rely on the results of Ref. [72] where 
a Boltzmann equation for the kinetic decoupling temperature T^d has been derived for general targets (but provided 
Mdm ^ m f)- From the Boltzmann equation given by Eq.(10) (and the results in Appendix A) of [72] we infer that 
departure from kinetic equilibrium between the scalar DM candidate and a specific fcrmionic species occured when a 
kinetic collision term, c(T), was of the order of the expansion rate H(T), 

~i/2 A 2 m 2 r 
M DM c(T kd ) ~ 9c f -J- \ dkk 5 u(k)e- u W r ™ ~ H(T kd ). (15) 



(2n)*M DM T kd 



,2 a „J =1/2 _ 1/2, 



Here A is the scalar-Higgs coupling, m/, the Higgs mass, u(k) — Jk 2 + m ? and g c ^ = g e g (T)/(l + T g' cS /Ag c ff) 

where g e e is the number of degrees of freedom contributing to the energy density at kinetic decoupling. 5 We find that 
kinetic decoupling occured at T kd ~ 300 MeV in the case of a 10 GeV candidate scattering with b/b quarks, while we 
find Tkd ~ 100 — 170 MeV for c/c quarks and r + /t~ respectively in the case of Mum — 5 and 10 GeV candidates, 
values which are not too far from our naive estimates. As for the muons, their coupling to the Higgs is parametrically 
smaller, but there are also lighter. In this case we find Tkd ~ 50 — 150 MeV, again for 5 and 10 GeV candidates. 
For lighter SM species, the decoupling occurs at much larger temperatures (we get for instance Ted ~ 700 MeV for 
e + /e~). Thus we take Tkd ^50 — 150 MeV as our estimate of the range of kinetic decoupling temperatures of the 
scalar DM candidate interacting through the Standard Model Higgs. 6 

In the case of "Z-like" we expect a priori Tkd to be ~ 1-10 MeV from previous estimates, like those based on the 
neutralino, since in this case the collisions with u/u and e + /e~ from the thermal bath are most important. To begin 
with we refer to the analysis of Ref. [71] , where an explicit formula for Tkd is given for the case of scattering of WIMPs 
on light, relativistic species in the thermal bath. Concretely, the thermal averaged elastic cross-section on relativistic 
light fermions is expanded as < er c i >= &q 1 (T/Mum) 1+1 and the dependence on the temperature is parameterized 
by I. In the framework of Ref. [71] kinetic decoupling is assumed to take place when the rate of collisions is of the 
order of the expansion rate, hence this approach is similar to the one we have advocated when deriving our naive 
estimates of the kinetic decoupling for the scalar interacting through the Higgs. In the Dirac fermion case, the elastic 
cross-section between the DM fermion and the SM fermions is mediated by the Z', is S-wave and proportional to 
dg 1 ~ sin <f> 2 G^M DM rnt^/m%,. The factor sin</> of suppression with respect to the elastic cross-section mediated by 
a Z boson is proportional to the kinetic mixing parameter e and to the mass of the Z' boson, as defined in Ref. 75j . 
Substituting CTq 1 in Eq.(4) of [74], the kinetic decoupling temperature is given by 

T kd ~ [^(^)M Pi ^sin^T 1/3 . (16) 

L -k z x Sir-* ' 1X1%, J 



4 



The Boltzmann equation for the kinetic decoupling temperature derived in |72| is valid for both relativistic and non-relativistic targets, 
but rests on an expansion of the collision term to leading order in parameters like uj/Mdm- One should keep in mind that this 
approximation may not be very good for the heaviest targets we consider, in particular bb quarks but, for the case of "Higgs-like" 
couplings, we only nee d to assume that it applies for collisions with lighter, albeit non-relativistic targets, which seems reasonable. 
The parameter g e g in jl5[ l may differ from g G ff by a factor of 0(2) between T = 100 and 200 MeV (see Figure 1 in Ref. 72 ), an effect 
we have however neglected in our order-of-magnitude estimates for T^d- 

We notice that these estimates imply that our candidates decoupled around the QCD phase transition. Beside the fact that one may 
be uncomfortable with considering free quarks below T ~ 1 GeV (an issue we do not address here), the QCD phase transition may 
also have an interesting impact on the formation of the smallest dark matter halos, and issue that has been addressed in, for instance, 
Ref. |73| . A relevant effect of the QCD phase transition is that there is a sharp drop of the speed of sound in the baryon-photon fluid. 
For dark matter that has kinetically decoupled before the QCD phase transition, this may lead to the formation of dark matter halos as 
light as 10 _1<) Mq or even less [73| . Whether this effect may be relevant for the candidates with "Higgs-like" couplings is unclear, given 
the proximity of the QCD phase transition critical temperature T c ~ 170 MeV and the kinetic decoupling temperature, so we leave this 
interesting possibility for a possible future study. 



Model 


T kd (MeV) 


M min 


(Mo) 


"H-likc" 


50-150 


5-10" 5 -3-10" 6 


8 • 10" 7 - 3 • 10" 8 


"Z-like" 


5- 25 


1- 10~ 3 -5- 10~ 5 


8- 10~ 4 -6- 10~ 6 



Table I: Range ol kinetic decoupling temperatures and the corresponding minimal halo masses, M m i n , for "Z-like" and "Higgs- 
like" DM candidates (see text). The extreme values for M m i n correspond to candidates with mass of 5 GeV (lower) and 10 
GeV (upper). Regarding M m i n , the results from Eqs. 1 17 1 and (181 are reported respectively in column 2 and 3. 



Notice that Tkd does not depend on the DM mass but depends on the details of the Z' model, specifically on the 
mass of the hidden gauge boson and the kinetic portal. In our estimation we use mz> = 15 GeV and e in the range 
5 x 10~ 4 — 5 x 10" 3 , values which have been shown in Ref. [75] to give the observed relic abundance for the light DM 
candidate. In this case T^d is 5 - 25 MeV, as reported in Table I, slightly colder than in the case of the singlet scalar, 
Tkd ~ 50—150 MeV. For the sake of comparison, we have also tentatively estimated the kinetic decoupling temperature 
using the equivalent of Eq. (15) for a Dirac fermion interacting through a Z' and found reasonable agreement, with 
21d ^5—15 MeV, hence we adopt here the range Tkd ^5 — 25 MeV for our estimates of the minimum halo masses 
for the Dirac fermion DM. 



From the kinetic decoupling temperature, one may estimate the free streaming length kf a and the mass of the lightest 
bound structure, or dark matter halo, Af min = 4/3 pdm(t /kfs) 3 Mq. We use here the results from Refs. [7H [71] and 
Ref. [72], which give slightly different answers. It is beyond our scope to choose between these options, so we give 
estimates for both cases. From the analysis in Ref. [72], one has 
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(17) 



Using this expression, we obtain M m ; n ~ 10 4 Mq for a "Z-like" candidate, and M m ; n ~ 10 5 Mq for "Higgs-like" 
couplings (see Table I). For an alternative approach, we consider the results of Ref. [76], and in particular 



M min «10- 4 f Tkd ) 
V 10 MeV/ 



M, 



o i 



(18) 



an expression which was obtained based on a numerical evaluation of kf s . Using our estimates for the kinetic decoupling 
temperature, we obtain halo masses that are typically lighter for "Higgs-like" couplings, possibly as small as Af min ~ 
10~ 8 Mq, than for "Z-like" couplings, with Af m j n ~ 10 -6 Mq respectively. We resume in Table I typical possible values 
of Ted and M m [ n for a Mom — 5 GeV and 10 GeV DM candidates, both with WMAP abundance. To emphasize 
the possibility that the scalar may potentially form lighter halos, in the evaluation of the diffuse gamma-ray signals 
we consider three different minimum mass values: M m ; n = 10~ 8 , 10~ 6 M© and 10 -4 Af© for the model independent 
and for the singlet case, while only two values for the "Z-like" DM candidate, Af m ; n = 10~ 6 Mq and 10 -4 Af©. For 
these models the minimum mass allowed by acoustic oscillations is far below the smallest halo mass permitted by 
the free streaming length scale fcf s [72]. Also, it is perhaps worth keeping in mind that some DM models designed 
to explain CoGeNT/DAMA may require different couplings than that of the Standard Model Higgs or the simple Z' 
model discussed here, opening the possibility of an even higher kinetic decoupling temperature, and thus lighter dark 
matter halos. 



We finally show in Figure [3] the boost, Eq. (11), and total enhancement factor, Eq. ( fl2| ), as a function of the 
redshift z. The (green) thin lines are for the Cwmap concentration parameter, while the (black) thick curves match 
the simple power law behavior for c V i r . The long-dashed lines correspond to M min = 10~ 8 Af©, the solid curves to 
•Afmin = 10~ 6 Mq and finally the short-dashed to Af m i n = 10 -4 Af© (the same convention for the curves description is 
maintained throughout the paper). On the left panel, we see that the enhancement due to structure formation starts 
to dominate over the smooth component at z < 30. The boost depends crucially on the concentration parameter, as 
Cwmap leads to smaller boosts with respect to Cpl; this has already been shown in [77], and we found agreement 
with their results. In the case of the Cwmap parameterization we have checked that our results are equivalent to 
those in Fig. (1) of Ref. [52]. On the right-hand side we show a zoom of the total dark matter density at small 
redshifts, the same quantity represented in Fig.(l) of Ref. [33]. The Fermi-LAT collaboration has used a different 
approach to compute the dark matter halo distribution, however we may compare their analytical estimation (labeled 



BulSub) with our calculations with M n 



10 Mq: for Cwmap we are 10% below their estimation, while for 
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Figure 4: Diffuse photon emission from a 10 GeV dark matter particle for BR b i — 100% on the left, and with annihilation 
into BR T + T - = 100% on the right. For the bb channel < av > a nn = 2.6 10 _26 cm 3 s _1 , while for the t + t~ case the thermal 
cross-section is fixed to 3.3 10 _26 cm 3 s _1 . The color code is as in Fig. pi The red points (below) are the measurements of the 
diffuse emission by Fermi-LAT [51] and the blue points (above) are those by EGRET J801I8T] . 



Cpl our boost is a factor ~ 3 above their curve. It is this small difference in the boost value at small redshifts 
(which are the values that contribute the most to the diffuse photon flux for a light WIMP) that will explain our 
more stringent upperbounds. We note that for a given concentration parameter and M m [ n = 1O~ 6 M0, the value 
-Mmin = 10~ 8 (10 -4 ) Mq corresponds to increase (decrease) the boost by a factor ~ 3. Even though the diffuse photon 
flux is also increased or decreased by the same amount, below we show the upper bounds for all the M m - ln cases and 
for the two fit of the c V i r . 

One should keep in mind the uncertainty related in the choice of profile, here a NFW. Less cuspy profiles, like 
Einasto [75] or Burkert [79] . will give a smaller boost, typically by a factor of 3 in the former case, and a factor of 8 
in the latter 1521. 



B. Discussion on the upper bounds 



The measured fluxes 0j GRB and the errors <Ji for the diffuse emission are those given by the Fermi-LAT collaboration, 
listed in Table I of [51 in the energy range from 0.2 GeV up to 102.4 GeV. The measurements arise from a full sky 
fit, for latitude b > 10° and from the data of the initial 10 months of the Fermi mission. From the isotropic diffuse 
component measured by Fermi, the IGRB flux is obtained subtracting the diffuse galactic emission, point sources 
and a cosmic rays background. The measured IGRB spectrum is compatible with a featureless power law with index 
—2.41 [ST]. In this light we derive conservative upper bounds on the annihilation cross-section of a WIMP, namely 
the dark matter signal should not exceed the measured flux in any individual bin i by more than a given amount 
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(19) 



The index n indicates the significance of the upper bound, n = 1.64 corresponds to the 95% confidence level (C.L.), 
assuming that the probability distributions for the intensity in each bin are Gaussian and independent from each 
other. The limits on the DM annihilation rate are derived taking the lowest value of < av > ann for which Eq. (19) is 
not satisfied. In order to explain the Fermi IGRB data one needs to add a background contribution in each energy 
bin where the DM flux is too small to account for the measured intensity. 

In Figure |4j the predicted IGRB photon flux times E 2 as a function of the photon energy is shown for a 10 GeV 
DM particle that annihilates with a branching ratio of 100% into bb or t + t~ together with the Fermi-LAT data (red 
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while it is fully compatible for M m ; n = 10 - ivi@ aim c vir — l^pl or ivi m - m iu ~iviq anu c„j r = u^map- 

shape of the flux, we see that the upper bounds will come from the first, second or at most third bin (depending 
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Figure 5: Fermi-LAT upper bounds in the plane < crv >„„„ versus the DM mass Mom at 95% C.L. for two different annihilation 
channels: on the left BR h i = 100% and on the right BR T + T - = 100%. The horizontal red (gray) region corresponds to the 
WMAP-7yrs 5 — a for the thermal annihilation cross-section. The vertical dotted line indicates the energy threshold for b, b 
and t production, on the left and right respectively. The color code is as in Fig. I3J 

on the dark matter mass) of the Fermi-LAT data, namely from measured photon energies lower than 2 GeV. The 
t + t~ case is analogous, again M min = 10~ 8 M Q and c„ r = Cpl excludes a cross-section of 3.3 • 10 _26 cm 3 s _1 , while 
this value is compatible for all other configurations of the astrophysical parameters. However in this case the most 
constraining energy bins are from the second to the fifth: the relevant energy range is more close to the photon 
production threshold, namely the dark matter mass. The contribution to the extra galactic photon flux comes from 
low redshifts. 

Scanning over the mass range from 1 GeV up to 20 GeV, we obtain the upper bounds in the plane < av > an n vs 
the dark matter mass Mom, Figure [5j The pure bb final state case is the left panel, while the right plot is given for 
an annihilation into t + t~ . Comparing these two model-independent cases, we note that the leptonic final state give 
more stringent upper bounds, for the same assumption on the astrophysical parameters, due to the different shape 
of the photon spectrum. In example, the < av > ann values that lead to a relic abundance in the WMAP range, are 
always compatible with the Fermi data for BR b i — 100%, while they are excluded at 95% C.L. for BR T + T - = 100% 
and Mum < 7 GeV, with M m ; n = 10 -6 M Q and Cpl. As already noted, the different choice in M m i n moves up or 
down the bounds by a factor of ~ 3. The bounds from the bb are compatible with those in |3" H 133 ] . The pure t + t~ 
has not yet been discussed for such low DM masses. 

IV. THE SINGLET SCALAR WITH HIGGS PORTAL AND A DIRAC FERMION WITH A Z' 



One may say that they are essentially two classes of models of dark matter. To the first class belong models within 
which dark matter is a by-product, the secondary effect of a deeper goal. Such are the neutralino and the axion, and, 
to some extent, DM candidates in theories with extra dimensions (for a review see |82j). The other class encompasses 
the many models which start with the DM problem, sometime in a minimalist way, sometime with a broader scope 
in mind. The first class may be nobler, but, with the advent of new data, the bottom-up approach to DM has gained 
some momentum. In particular, the interpretation of the CoGeNT or DAMA data as being due to DM has motivated 
much works, both on the possibility to explain the data with a light neutralino [44] I83M88] and on new models (see 
e.g. [51 ItH 131)1 1751 l85TI5ni89"rl9T>] ). In this section we consider two very simple light WIMP candidates. These models 
have in common the facts that they are based on minimal extensions of the SM and yet may explain the direct 
detection data while having a relic abundance that agrees with cosmological observations. The first model is based 
on a real scalar singlet, which is interacting through the Higgs portal with the SM degrees of freedom [37, 100, 101 . 



The second model invokes a singlet Dirac fermion, 
gauge boson [73 [Ml ESI ■ 



interacting through the kinetic mixing portal, i.e. via a new Z' 
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Figure 6: Fermi-LAT exclusion limits (colours as in Fig.[5| on the singlet dark matter candidate with "Higgs-like" annihilation 
channels, together with the 2 — a (purple, dark grey) and 5 — a (pink, light grey) WMAP-7yrs regions. 



That a real scalar singlet interacting through the Higgs may be in agreement with both CoGeNT-DAMA and 
WMAP has been argued in [91 1571 HTH] . This is implicit in other works. In particular in the results of 8J, which 
are based on an effective operator approach (see in particular Figure 4, lower left plot). We may also refer to |98j . 
which is based on a singlet complex scalar, in particular Figure 7. In practice, these models differ only in the choice 
of couplings to the SM degrees of freedom, which in the case of [HI I37J ar e simply proportional to the SM fermions 
mass ("Higgs-like" couplings) and thus more constrained. These scalar DM models also have in common the fact 
that, at the low energies relevant for the CoGeNT-DAMA regions and the calculation of the relic abundance from 
freeze-out, the ratio of the annihilation and elastic scattering on nuclei cross-sections only depends on the DM mass. 
In particular, for the case of one real scalar singlet S of mass ms interacting through the Higgs, the ratio is 
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where n c — 3(1) for quarks (leptons), fi n is the DM-nucleon reduced mass and the factor / parameterizes the Higgs 
to nucleons coupling from the trace anomaly, which has central value / ~ 0.3 [37] . Taking into account the branching 
ratio for annihilation of a light S into SM fermions (see Table II), the limits on the annihilation cross-section of the S 
from the IGRB measured by Fermi-LAT are shown in Figure [6] Which limiting curve we may believe is most relevant 
depends strongly on the minimum halo mass M m i n . An interesting feature is that, after freeze-out, a particle with 
"Higgs-like" couplings to SM degrees of freedom, has very little interactions with particles from the thermal bath. 
Indeed we may neglect the Yukawa co upli ng with the electron and positrons, while the protons and neutrons are very 
few in the thermal bath. Using Eq. (18 1 we estimated M m - m ~ 10~ 8 M@. To this parameter, together with Cpl, 
corresponds the strongest constraint we may get based on the IGRB (long dashed curves in Figure [6]) . In this plot 
we see that all the < av > a „„ values that account for the WMAP abundance are excluded at 95% C.L., while more 
relaxed constraints come from either M min = 10~ 6 M Q and Cpl (solid black thick curve) or M min = 1O _8 M and 
Cwmap (long-dashed thin green line). Annihilation cross-sections of the order 3 • 10 _26 cm 3 s _1 are excluded only for 
very light scalar masses, ms < 7 GeV. 

Since the lowest energy bin measured by Fermi-LAT is most relevant for light WIMPs, one may ask whether the 
EGRET data, which extend to lower energies could give further constraints [HUllBTj . However, this is not the case for 
the range of mass we consider, as is shown in Figure [7] for a ms = 2 GeV candidate. A similar conclusion holds for 
the Dirac fermion candidate. Of course, this does not mean that EGRET data are not relevant for even lighter DM 
candidates, Mdm % 1 GeV, but only that such candidates are beyond the scope of the models discussed here. 



Using the Eq. ( 20 ) , we may transpose the limits on the annihilation cross-section into limits on the elastic scattering 
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Branching ratios 


m s 


bb 


cc 


T T 


others 


20 GeV 


85% 


5% 


9% 


~ 1% 


lOGeV 


83% 


7% 


10% 


< 1% 


5 GeV 


16% 


36% 


42% 


-5% 


2 GeV 


\ 


69 % 


22% 


-9% 



Table II: Branching ratios in the main annihilation channels of the scalar singlet, for various light candidates. 







E (GeV) 



Figure 7: Diffuse photon emission for a ms — 2 GeV candidate, the lightest DM particle we consider here, compared to 
the Fermi-LAT (in red, below) and EGRET data (in blue, above). Because the EGRET data points are above those from 
Fermi-LAT, only the latter is constraining for the mass range relevant for the model. From Figure |6| we used the maximum 



allowed < av > ann — 3.2 • 10 26 cm 3 s * and the color code as in Figure HI 



cross-section on nuclei, and plot them in the c° — ms plane, Figure p| on the left, together with the CoGeNT-DAMA 
regions and the exclusion limits set by CDMS-Si and XenonlOO, which were computed as in 9j- Following the same 
convention as in the other figures, the upper bounds from the IGRB are shown in green or black. For M m { n — 10 _8 M Q 
and a power law concentration parameter the DAMA and CoGeNT regions are both excluded at 95% C.L., while 
the same minimum halo mass with C^map gives limits which are marginally compatible with the CoGeNT region. 
A choice of M m i n = 10~ 6 M Q and Cwmap is totally compatible with the CoGeNT region and can be accommodated 
with the DAMA region (see footnote 7). 

For comparison, we also give in Figured (right panel) the region of the er"— ms plane consistent with WMAP. In this 
figure, the black continuous lines is based on the hypothesis that the QCD phase transition occurred at a temperature 
T c = 150 MeV (which is the standard assumption both for DarkSusy [105j and Micromegas |106| ). Freeze-out after 
the QCD phase transition would require an annihilation cross-section which is about twice that for a standard WIMP, 
< av >ann ~ 3- 10~ 26 cm 3 s _1 (see for instance [83]). Since typically Xf = Mdm/Tjo ~ 20, this effect is only relevant 
for Mjjm < 3 GeV. However, we may contemplate the possibility that the QCD phase transition took place at a 
higher temperature, for instance 150 MeV < T c < 500 MeV, which is relevant for a candidate Mum ^S 10 GeV, and 
thus for the CoGeNT-DAMA regions. 7 



7 This effect may be relevant in the light of the many possible uncertainties that may be hidden in a Figure like Isl (plot on the right). For 
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Figure 8: The SI cross-section (<r°) vs the scalar singlet mass (ms). On the left, the region (in green) corresponds to 
CoGeNT [3] (minimum % 2 , with contours at 90 and 99.9% C.L.). The black thick and green thin lines are as in Fig. pjl and 
denote the exclusion limits from the Fermi-LAT diffuse gamma-rays flux at 95% C.L.. On the right, we give, for the sake of 
comparison, the corresponding exclusion limits from direct detection experiments. We include the DAMA regions .2 (goodness- 
of-fit, also at 90 and 99.9% C.L.) both with (below CoGeNT, purple/orange) and without (above CoGeNT, purple, no fill) 
channelling. The blue (short-dashed) line is the 90% C.L. exclusion limit from CDMS-Si [15]. The red solid curve is the 90% 
exclusion limit from XenonlOO |19| . using LeffMin and a threshold at 3 PhotoElectrons. For more details on the experimental 
upper bounds and conventions see [9]. On the right, we include the envelope of a 1 ^ — ms parameters consistent with WMAP. 
The continuous region corresponds to the standard assumption of a QCD phase transition at T c = 150 MeV. The black dashed 
lines extend this domain for other, possible but less likely, values of T c , from left to right T c = 300 MeV and T c = 500 MeV 
respectively. 

One may notice that the constraints from IGRB extend to candidates which are even (slightly) lighter than those of 
interest for CoGeNT-DAMA, and which are not constrained by direct detection experiments (so far). In particular, 
given the conservative constraint corresponding to the continuous black curve in the left panel of Figure [51 we may 
conclude that candidates lighter than 5 GeV, and consistent with WMAP, are excluded at 95 % C.L. by the Fermi-LAT 
data on the IGRB. 

For the sake of comparison with the singlet scalar candidate, we consider a model that consists of a singlet Dirac 
fermion candidate ip, charged under a broken U(l)' gauge group, which may interact with the SM degrees of freedom 
through the kinetic mixing portal. In |75j . it has been shown that such a candidate may be consistent with both 
WMAP and the CoGeNT-DAMA regions as well as constraints from LEPI on the Z invisible width. We notice that 
this result is a priori in contradiction with the one drawn in [37] and [5] , where it has been shown that a singlet Dirac 
DM candidate (for a single, vector or scalar annihilation channel) that fits the CoGeNT-DAMA regions has a too 
large relic abundance. However, there is no magic, as the conclusion in [5] I37| has been reached by assuming that 
there exists a one-to-one correspondence between the annihilation and scattering cross-section (for fixed dark matter 
mass), while in [75] the proximity of the Z' pole is used to enhance this annihilation cross-section so as to get the 
right relic abundance. Incidentally, this will prevent us from presenting the exclusion limits in the a^ — Mdm plane 



instance, in 23 , it is argued that channelling may, after all, not be very relevant for the interpretation of the DAMA data, so the lower 
of the two regions corresponding to DAMA most presumably does not exists. There is however still substantial freedom just on the 
experimental side, and, based on the current uncertainties on parameters like quenching, the authors of 1141 argued quite convincingly 
that, not only, both the DAMA (without channelling) and CoGeNT may be consistent with each others but, moreover, are not excluded 
by the current exclusion limits. 
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Figure 9: The same as in Fig. [6] for the singlet Dirac fermion candidate, through annihilation channels with "Z-like" couplings. 



for this model, for the ratio of < <rv > an n and cr° depends on mz' and Tz> 
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In Eq. (21), af and vj are the SM model axial and vector couplings of the fermions to the Z. The limits on the 
annihilation cross-section of tp from the IGRB are shown in Figure [9| That the Z' has "Z-like" couplings to the SM 
degrees of freedom (compared to the "Higgs-like" couplings of the scalar singlet) has two interesting consequences 
regarding the constraints from the IGRB. First the annihilation channels are distinct, with a spectrum that is harder 
with respect to the scalar singlet one, and lead to the constraints from Fermi-LAT to extend to rather light DM 
candidates. Also, because of the abundance of the e + /e~ (and neutrinos) in the thermal bath, the kinetic decoupling 
(compared to the chemical freeze-out that fixes the relic abundance) of a DM candidate with "Z-like" couplings occurs 
at a much smaller temperature, which, following [73], we estimated to be O(10 MeV). This implies that the minimum 
halo mass that may form in t he mo del is substantially larger for the Dirac fermion than in the case of the singlet scalar, 
•^min ~ lO _4 A/ from Eqs^ jl7|l8 ). This is the reason why we only plot the exclusion limits for M m ; n ~ 1O~ 6 -M0 and 
M m i n ~ 1O _4 M in Figure|9j In this case the most stringent limit, given by M min = lO~ 6 A/ and Cpl excludes at 95% 
C.L. the WMAP region for m^ < 5 GeV. All the other curves only marginally affect < av > ann ~ 3 • 10 _26 cm 3 s _1 . 

Altogether, we may conclude that, although they are both viable candidates with respect to WMAP and CoGeNT- 
DAMA, the constraints from the IGRB are substantially weaker for the -0 candidate, with "Z-like" couplings, than 
for the S candidate with "Higgs-like" couplings, a feature we believe should be of interest for model building. 

Needless to say, such constraints should be taken with a grain of salt, as they suffer from the many astrophysical 
uncertainties, regarding for instance the profile of dark matter in halos, and other parameters which enter in the 
estimation of the extragalactic boost (even the use of the Press-Schechter formula may be a matter of discussion, as 
is the extrapolation of over many decades of the known power spectrum of the DM inhomogeneities) . Regarding the 
profile we have chosen use the somewhat standard NFW choice, like in [3H |33]; clearly a different choice, like the 
Einasto or Burkert profiles would give less stringent constraints. 



CONCLUSIONS AND PROSPECTS 



In this work, we have studied further the constraints that may be set, using the Fermi-LAT data on the isotropic 
diffuse gamma-ray emission, on the parameter space of two simple, albeit generic, DM candidates, which has been 
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shown to be consistent both with the WMAP relic abundance, and the regions favoured by CoGeNT-DAMA [51 175] . 
The constraints we get are consistent with results which may be inferred from other current works on constraints 
on DM from extragalactic gamma-ray fluxes [3TJ [33J . We have however extended them to lower DM masses and 
make explicitly the connection with the CoGeNT-DAMA results tentatively interpreted as being due to the elastic 
scattering of DM particles. Our main results are summarized in Figure [6] for the case of a singlet scalar candidate 
S interacting through the Higgs portal ( "Higgs-like" couplings to SM degrees of freedom) and Figure [9] for the case 
of Dirac fermion candidate ij> interacting through a Z' ("Z-like" couplings). In the former case, we have used the 
one-to-one correspondence between the elastic cross-section and the annihilation cross-section, Eq. ( |20[ ) to express the 
limit from the isotropic gamma-ray background radiation (IGRB) directly in the cr° — ms plane, Figure pi which is 
the most direct way to express the constraints on CoGeNT-DAMA from the IGRB. This works only for the scalar 
candidate, because it has so few free parameters, but would also apply to an approach based on effective operators. 

Like all indirect signatures, the constraints on DM derived from the data on the IGRB are quite sensitive to 
astrophysical uncertainties. As such, they are perhaps best thought as being complementary to those that could be 
extracted from other observations. Most relevant are the constraints based on gamma-rays from dwarf spheroidal 
galaxies (dSphs), also observed by Fermi-LAT, which may also exclude the light WIMPs, Mjjm ^ 10 GeV, assuming 
a NFW profile of DM distribution [8, 91 1107] . Unlike the constraints from the IGRB, the constraints from dSphs 
are the same for the scalar singlet and the Dirac singlet, as they both have S-wave, unsuppressed, annihilation into 
SM particles. A potentially important indirect constraints on a light WIMP could come from measurements of the 
synchrotron radiation near the galactic centre, and in particular from galaxy clusters ;46^50j . Such constraints are a 
priori the strongest for candidates that annihilate directly into e + e~ pairs (see e.g. [50]), and so we expect that they 
may be stronger for the Dirac singlet interacting through a Z' than for the singlet scalar candidate discussed here. A 
multi- wavelengths analysis of the light WIMPs mass range could thus be of interest, but this is beyond the scope of 
the present work. 
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